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PACS 82 . 70 . Dd - Colloids 
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PACS 75 . 75 . Jn - Dynamics of magnetic nanoparticles 

Abstract - In colloidal suspensions, self-organization processes can be easily fueled by external 
fields. One particularly interesting class of phenomena occurs in monolayers of dipolar particles 
that are driven by rotating external fields. Here we report results from a computer simulation 
study of such systems focusing on the clustering behavior also observed in recent experiments. 
The key result of this paper is a novel interpretation of this pattern formation phenomenon: We 
show the clustering to be a by-product of a vapor-liquid first order phase transition. In fact, the 
observed dynamic coarsening process corresponds to the spindodal demixing that occurs during 
such a transition. 
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Introduction. — Self-assembly and self-organization 
rocesses of colloidal particles are topics that have been 
receiving much attention in the recent past. 

Indeed, such systems display a multitude of equilibrium 
and non-equilibrium self-assembled structures, examples 
being lane formation [TJEJ, shear banding [3J, the coiling 
up of magnetic chains [][] , and the wide range of patterns 
observed in particles immersed in liquid crystals [5] . Here, 
we are particularly interested in colloidal systems involv- 
ing dipolar interactions. Prime examples of the resulting 
self-assembled structures include chain formation in con- 
stant external fields [BHH] , layer formation in rotating fields 
[7HTU] , and the structure formation of colloidal particles in 
triaxial fields [HHI3]- 

Up until recently, most works on self-organization un- 
der the influence of time-dependent external fields focused 
on induced dipolar particles. A noteworthy exception is a 
paper by Murashov and Patey, in which layer formation 
of rotationally driven colloidal particles carrying a perma- 
nent dipole moment was investigated [14) . 

A particularly interesting self-organization process oc- 
curs when colloidal particles are exposed to rotating fields 
in a quasi-two-dimensional geometry. In this situation, 
the external fields arc found to induce the formation of 
two-dimensional clustered structures. Not only does this 
work for particles in which a dipole moment can be in- 
duced [T5HT7] , but also for particles carrying a permanent 
dipole moment. This was recently shown by Weddemann 



and coworkers in an experimental study (T5II19J . 

In the present paper, we want to pick up on this phe- 
nomenon and provide a novel interpretation of the two- 
dimensional cluster formation process. Specifically, we will 
show that this self-organization process is a by-product of 
an equilibrium liquid- vapor phase transition. 

We investigate the colloidal system by making use of 
different computer simulation techniques. In these simu- 
lations, we use dipolar particles in a quasi-two-dimensional 
geometry to model the system. This means that the 
dipoles can rotate freely in all the spatial directions while 
the translational motion is restricted to a two-dimensional 
plane. We use Brownian dynamics simulations to under- 
stand the dynamical properties of the system and Wang- 
Landau-Monte-Carlo simulations to look into its phase be- 
havior. Additionally, to assess the influence of the solvent 
on the system, we take it implicitly into account by em- 
ploying Brownian dynamics simulations that include hy- 
drodynamic interactions. 

This paper is organized as follows: After introducing the 
model and the different simulation techniques, we first dis- 
cuss the full non-equilibrium "phase" diagram indicating 
the region of cluster formation in the domain of frequency 
and strength of the external field at constant equilibrium 
thermodynamic parameters. In a next step, we investigate 
the influence of hydrodynamic interactions on the forma- 
tion of clusters. Then we present the principal point of 
this paper: We show that the non-equilibrium cluster for- 
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mation is essentially an equilibrium phase transition. To 
do this, we calculate an equilibrium phase diagram, in 
the construction of which an effective non-time-dependent 
inter-particle interaction is used, and examine the growth 
of the characteristic domain size of the clusters. The paper 
is then closed with a brief summary and conclusions. 

Model. — To model the (dipolar) colloidal particles 
in our simulations we use a dipolar soft sphere (DSS) po- 
tential, which is comprised of a repulsive part U Tep and a 
point dipolc-dipolc interaction part 



ff Ubb (r y , /«<,/*,) 
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In eq. (HJ) , Yij is the vector between the positions of the 
particles i and j, r^ its absolute value, and \x i is the dipole 
moment of the ith particle. The potential U rcp is the 
(shifted, cf. [50]) soft sphere (SS) potential for particles of 
diameter er, with the unshifted SS potential being given 
byU ss (r)=4e(a/ nj ) 12 . 

We investigate the system by making use of different 
simulation techniques. First, we employ non-overdamped 
Brownian dynamics (BD) simulations. The corresponding 
equations of motion for particles of mass m and moment 
of inertia I are [T4"ll2"l"] 
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where £t and £r are the translational and rotational 
friction coefficients, respectively. Consistent with ear- 
lier simulation studies of fcrrofluidic particles [TU1IH] . 
£t = 13.5-^/me/cr 2 and £r = 0A5Vrnea 2 were used. Fur- 
thermore, Ui is the angular velocity of particle i, Ff and 
Tp are random Gaussian forces and torques whose vari- 
ance is related to the friction coefficients [14] , and T° xt are 
torques due to an external field. The field is homogeneous, 
rotates in the dipolar monolayer, and is given by 



B(i) = Bo(e x cosuot + e y sinwoi), 



(4) 



where luq is the frequency of the field and Bq its strength. 

The equations of motion @ and (J3) were integrated 
with a Leapfrog algorithm [50] using a time-step of At = 
0.0025(mcr 2 / £ ) 1/2 and 4900 or 1225 particles, respectively. 

Further, to investigate the influence of a solvent within 
our implicit model, we use a Brownian dynamics simu- 
lation that includes hydrodynamic interactions between 
the particles. These interactions are incorporated up to 
third order in the diffusion tensor for the translation- 
translation coupling, the rotation-rotation coupling, and 
the translation- rotation (and vice versa) coupling [22j[23] . 
The tensor describing the translation-translation coupling 
is the well known Rotne-Prager tensor. The other cou- 
plings incorporate the influence of the additional rota- 
tional degrees of freedom. The time evolution equations 



that were used can be found in |23| and we used 324 par- 
ticles in these simulations. 

Finally, to investigate the equilibrium phase behavior 
of the system (based on a time-averaged potential, see 
cq. (O), we use Monte Carlo (MC) simulations in the 
Grand-canonical ensemble. In general, first-order phase 
transitions arc plagued by a large free energy barrier sep- 
arating both phases, making unbiased sampling very inef- 
ficient. In order to overcome the barrier we use a method 
proposed by Wang and Landau |24j which significantly 
improves the sampling by adding a weight-function to the 
simulation. Details of the method applied to dipolar sys- 
tems in two dimensions can be found in |25j . 

In the case of BD and MC simulations we deal with the 
long-ranged dipolar interactions by using the Ewald sum- 
mation method [55]. In our hydrodynamic simulations, on 
the other hand, we only consider a single simulation box 
filled with dipolar particles. 




Fig. 1: Snapshots showing a system in zero field (a) and ex- 
posed to a rotating field of strength (e/cr 3 ) 1 ' 2 ^ = 50.0 and 
frequency (ma 2 /e) _1 ' 2 cjo = 20 (b). Snapshot (b) was taken at 
time t — 375(mcr 2 /e) after the start of the simulation. 4900 
particles were used. 

Synchronization and cluster formation. A 

snapshot of a quasi-two-dimensional system (pa 2 = 0.3, 
ksT /e = 1.0) of strongly coupled dipolar particles (Add = 
^/ksTcr 3 = 9) at zero field is shown in fig.QJi. As is typ- 
ical for such dipolar systems the particles align in a head- 
to-tail configuration, which, in a two-dimensional geome- 
try, results in the formation of chains and rings |27j . 

If we expose such a system to a rotating in-plane field of 
sufficient strength and frequency, we observe the particles 
agglomerate into two dimensional clusters. An example of 
this can be seen in fig.[T}D. The observed clustering behav- 
ior already indicates that there are attractive interactions 
in the system, which play a crucial role. The origin of 
these interactions can be rationalized as follows: Averag- 
ing the dipolar inter-particle interaction potential over one 
rotational period of the field under the assumptions that 
the particles do not move translationally and rotate syn- 
chronously with the field (i.e. follow the field at constant 
phase difference) yields 



u \ r 'j> - 9? ,3 • 



(5) 
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Clearly, for (|5|) to be a good approximation to the true 
inter-particle interaction, it is crucial that essentially all 
the particles follow the field. An extensive analysis of this 
synchronization behavior of the dipolar particles with the 
field in three dimensions can be found in [TCJIl2"5] . 

To systematically investigate the appearance of synchro- 
nization and clustering we scanned a wide range of fre- 
quencies and field strengths. We consider a system as 
clustered if the particles have on average more than 2.3 
neighbors within a distance of 1.7a from their center. The 
latter value was used, since it is slightly larger than the 
typical distance between neighboring particles in a clus- 
tered system. This was found by looking at the first 
minimum of respective pair correlation functions. The 
fact that we require 2.3 neighbors on average ensures that 
two-dimensional aggregates are counted as clusters while 
chainlike structures are disregarded. 

The results of this investigation of the space of the 
field parameters can be seen in fig. [5J Depicted are 
three distinct regions, denoted "synchronous" , "syn- 
chronous/clustered" , and "not synchronous" . The first 
region is comprised of systems in which the particles ro- 
tate synchronously with the field but do not form clus- 
ters. Within this region, chains in the direction of the 
field can be observed at low frequencies while spatial in- 
homogeneities begin to appear at larger frequencies. As 
also becomes apparent here, a minimal Bq is required for 
the field to align the particles with itself. 

In the second region, the synchronous rotation con- 
tinues but is now accompanied by the formation of 
two-dimensional clusters. This indicates that the effec- 
tive inter-particle potential becomes sufficiently isotropic 
within this region to be reasonably described by the aver- 
aged dipolar potential ([5]). 

As expected, the diagram in fig. [5] also depends on the 
friction coefficients used in the BD simulations. In particu- 
lar, we found the size of the synchronous regime in the fre- 
quency domain to shrink in favor of the "clustered" regime 
if the translational friction coefficient £t is increased. This 
results from reduced displacement of the particles after one 
rotational period of the field thereby making eq. ([5]) valid 
at much lower driving frequencies. 

In the third region, we find neither synchronization nor 
cluster formation. Clearly, the lack of synchronization is 
the direct cause of the breakdown of cluster formation 
(cf. eq. ©). The loss of synchronization occurs, since the 
torques on the particles due to the external field become 
unable to overcome the torques due to the rotational fric- 
tion (similarly to what is seen in three dimensions |10jb 

Given the clustering scenario and our explanations so 
far, it is important to ask to which extent the BD sim- 
ulations can describe the dynamics of the real colloidal 
system. The latter includes, by definition, a solvent. The 
rotating external field, which constantly generates rota- 
tional motion of the particles, will create flow fields that 
can result in considerable motion of the particles. This 
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Fig. 2: Synchronization behavior and cluster formation de- 
pending on frequency ojo and strength Bo of the field. The 
density and temperature of the system used were pa 2 = 0.2 
and ksT/e — 1.0, respectively. 1225 particles were used is 
these simulations. 



might influence the cluster formation phenomenon due to 
the following reason: The averaged potential ([5]) is only 
valid as an approximation to the true inter-particle inter- 
action if the translational motion of the particles during 
one rotational period of the field is small. 

Hence, in order to find out whether cluster formation 
persists when hydrodynamic interactions are present, we 
built a simulation that takes these into account. To assess 
their influence on the cluster formation, it is sufficient to 
consider only a single simulation box filled with particles. 

We considered a number of state points (pa 2 = 
0.1,0.2,0.3, k B T/e = 1.0,1.8 at ua 2 /D = 250, 
Bo(a 3 /e) 1 ^ 2 = 80) in our simulations with and without 
hydrodynamic interactions included. We found that clus- 
ter formation occurs in both these cases at all the con- 
sidered densities and temperatures, with a single cluster 
eventually forming in the simulation box. This can be 
seen in fig. [31 where we show snapshots of the evolution 
of a rotationally driven system. The snapshots in the top 
row show a system in which hydrodynamic interactions are 
not taken into account, while hydrodynamic interactions 
are present in the snapshots in the lower row. Another 
important point that is illustrated by fig. [3J is that the 
cluster formation process is considerably accelerated by 
the hydrodynamic interactions. At the intermediate time 
(£' = 10.8), only a single cluster remains in the hydrody- 
namically interacting system (fig-EJ;), while it takes much 
longer for the not hydrodynamically interacting system to 
reach the same state (cf. fig. [3b, c). 

Accelerated cluster formation is a direct consequence of 
both the translation-translation and rotation-translation 
coupling. In particular the former coupling accelerates the 
formation of clusters. This can be inferred from selectively 
switching off the different hydrodynamic couplings. Fur- 
ther, note that if hydrodynamic interactions are present, 
the cluster rotates rapidly around its center, which is not 
the case if these interactions are absent. Here, the sole 
cause is the hydrodynamic rotation-translation coupling. 
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Fig. 3: Snapshots of a system at different times without (top) 
and with (bottom) hydrodynamic interactions. The snapshots 
in the first column were taken at t' = 2.5, the ones in the 
second column at t' = f0.8, and the one in the last column 
at t' — 1 38.5 after the start of the simulation. Note that t' = 
tDo/a 2 



Thus, the hydrodynamic interactions affect the (rota- 
tional) cluster motion but do not hinder the particles in 
the cluster formation process at all. 

Relation to condensation. Given the clustering 
behavior of the non-equilibrium, yet fully synchronized, 
field-driven system, we now ask to which extent the be- 
havior of the system can be understood by that of an 
equilibrium system interacting via the effective interac- 
tion given by cq. (|5|) (as well as the repulsive potential 
rjrep ^ c £ C q^ (PQ))). More precisely, our hypothesis is that 
the observed cluster formation stems from an equilibrium, 
first-order phase transition between a vapor and a liquid 
phase. To test this hypothesis, we have performed Wang- 
Landau-Monte-Carlo simulations of a system interacting 
via U ID (eq. ©) [3QJ. 
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Fig. 4: Phase diagram in the T-p domain at (ecr 3 ) 1 ' 2 ^ = 3. T c 
denotes the critical temperature, the solid line corresponds to 
the binodal representing coexisting gas and liquid states. 

The key result we obtained from these simulations is 
the phase diagram presented in fig. 0] The presence of 
a binodal shows that there is indeed a gas-liquid phase 



transition in the non-driven system. 

Moreover, if we now compare this phase diagram 
with the thermodynamic state point related to our non- 
equilibrium system (see e.g. fig. [T})), we find that the lat- 
ter {pa 2 = 0.3, ksT/e = 1) lies well within the two-phase 
region. This finding clearly suggests that the cluster for- 
mation observed in the (fully synchronized) field-driven 
system is spinodal decomposition, i.e., a dynamic phe- 
nomenon occurring during a first-order phase transition. 
To check this conjecture in more detail, we looked at the 
time evolution of the cluster sizes £. For phase separating 
systems it is well established that £ exhibits power law be- 
havior |31j . i.e. £ oc t a , with the corresponding exponents 
depending on the growth stage. Such a behavior is also 
seen in Molecular dynamics (MD) simulations. In particu- 
lar, domains with growth proportional to t 1//2 J33JI33] and 
t [Mj have been identified. 

These power laws are universal in MD simulations but 
they do not necessarily apply to BD simulations with their 
modified equations of motion. This was shown, e.g., by 
Lodge and Heyes for the case of overdamped BD [3S] . At 
the same time, however, the clusters in ref. [3S| were still 
found to grow with a power law. To check for the existence 
of cluster growth with a power law in non-overdamped BD 
simulations, we first investigated a "reference system", 
whose equilibrium behavior is well studied. Specifically, 
we considered a two-dimensional Lennard- Jones system 
at a 2 = 0.3, ksT/e — 0.45 with the critical point being 
at pa 2 w 0.335, k B T/e w 0.533 [M]. Investigating the 
domain size, we did indeed find a power law dependence 
£ oc t a with a w 0.30. Note that the cluster size was 
obtained by measuring the distance at which the pair cor- 
relation function assumes a value of one for the first time 
if the radial bins are taken to be larger than the particle 
diameters (cf. [37]). 

Similarly, we checked the cluster growth for a driven 
dipolar system and a system interacting via the potential 
Z7 ID . In the simulations we used a density and a temper- 
ature of pa 2 = 0.3 and ksT/e = 1.5, respectively, which 
put the systems well inside the coexistence region of fig. |U 
The domain sizes over time that were extracted from the 
simulations are shown in fig. [5j As can be seen, the clus- 
ter sizes of these two systems are almost identical at any 
given time. In particular, the characteristic domain sizes 
grow with a power law £ oc t a with a being equal to 0.36 
in both cases. We note that this is almost identical to the 
Lifshitz-Slyozov growth law {£ oc i 1 / 3 ) [31j . 

From these two results we conclude that the non- 
equilibrium system does indeed undergo spinodal decom- 
position. First, the cluster growth proceeds with a power 
law, which is typical within the spinodal region. Second, 
the growth behavior remains unchanged even if the in- 
teractions between the driven dipoles are replaced with 
the effective ones. This emphasizes the similarity between 
those two systems and once more shows the validity of the 
phase diagram for the driven system. 

Note that hydrodynamic results are not included in 
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Fig. 5: Cluster growths for a system interacting via the effective 
potential and a system that is driven by an external rotating 
field at density pa 2 = 0.3 and temperature fcsT '/e = 1.5 (state 
inside the coexistence region of fig. [4]). 4900 particles were used 
in the simulations. 



fig. [5] We only considered 324 particles in these simu- 
lations and additionally did not use periodic boundary 
conditions. This makes an accurate determination of the 
cluster growth very difficult. 

Conclusions. — In this paper we have investigated 
the formation of two-dimensional aggregates in monolay- 
ers of dipolar particles that are driven by rotating external 
in-plane fields. 

The first result of this paper is a non-equilibrium 
"phase" diagram, which shows the regions of synchroniza- 
tion and cluster formation in the u — -Bn-domain. At high 
frequencies the synchronization of the particles with the 
field breaks down, which results in a breakdown of clus- 
ter formation. Similarly, we do not find cluster formation 
at low frequencies. Here, the effective interaction between 
the particles is not well enough described by C/ ID , since the 
particles move considerably during one rotational period 
of the field. This changes in-between those frequencies, 
where the particles rotate synchronously and sufficiently 
fast, which leads to the formation of clusters. 

Next, we investigated the stability of the clustering phe- 
nomenon when hydrodynamic interactions are present. In 
our simulations, we found the phenomenon to persist de- 
spite these interactions. In fact, the cluster formation 
seems to proceed at even faster rates with these inter- 
actions included. We attribute this to a combination 
of the hydrodynamic translation-translation and rotation- 
translation coupling. The former seems to accelerate the 
formation of clusters while the latter leads to the quick for- 
mation of a single cluster in the center of the simulation 
box. 

Test simulations indicate that these interactions 
have additional consequences: Compared to a non- 
hydrodynamically interacting system, they seem to allow 
for the formation of clusters at lower driving frequencies 



and therefore affect the non-equilibrium "phase" diagram 
in fig. [5J To study the precise influence of the hydrody- 
namic interactions is outside of the scope of this paper, 
but would be very illuminating in its own right. 

We concluded our analysis with the main result of this 
letter: We established the clustering phenomenon to be a 
by-product of a liquid-vapor phase transition. This was 
done in two steps: We began by uncovering a phase tran- 
sition via Wang-Landau-MC simulations in a system in- 
teracting via the effective potential U lB . Recall that this 
potential describes the interactions between the particles 
very well in the driven system at sufficiently high frequen- 
cies. In a next step we examined the domain growth of 
the driven system within the binodal region of the phase 
diagram. As expected within the spinodal region of a 
liquid-vapor phase transition, we found the characteris- 
tic cluster size to grow with a power law. Additionally, 
it essentially agrees with the domain growth of the non- 
driven system interacting via the effective potential [/ ID . 
These facts lead us to conclude that the clustering process 
corresponds to the pattern formation occurring inside the 
binodal of a vapor-liquid phase transition. 

Given these findings, it is interesting to compare them 
to recent experimental results. Indeed, cluster formation 
in monolayers resulting from a rotating external field has 
been observed multiple times [15HI9] . In most of these 
publications induced dipolar particles are brought to self- 
assemble into two-dimensional aggregates. The only paper 
in which particles with a permanent dipole moment were 
used (ref. [IB]) features a dipole-dipole coupling strength 
(Add = l^/kBTa 3 ) that is dominated by the dipole-ficld 
coupling strength (Adf = ^Bo/hsT). There, the ratio 
Adf/Add is about 6, which is larger than the largest ratio 
that appears in fig. [2] Consequently, we expect the parti- 
cles in [T5J to rotate synchronously with the field, resulting 
in the effective interaction [/ ID and the observed cluster 
formation. 

The clusters found in the literature are typically hexag- 
onally ordered. In our simulations this becomes more and 
more true with increasing frequency and strength of the 
field as well as with decreasing temperature. It should 
be noted, however, that our BD simulations involve large 
Brownian kicks, which result in deviations from the hexag- 
onal structure. In the references P^llTTrH^] . micrometer- 
sized particles were used, which makes these contributions 
much less significant leading to more pronounced order. 

Further, the driving frequencies used in these publica- 
tions are considerably smaller than the ones used here. 
This can once again be explained by the size of the parti- 
cles: Larger particles typically have smaller friction coef- 
ficients, which, as test simulations show, result in cluster 
formation at lower frequencies of the field. 

The coarsening process investigated in this paper can 
universally be observed in phase separating systems. Con- 
densation transitions and spinodal demixing in binary flu- 
ids or metallic alloys are popular examples of this. The 
process reported in this letter is exceptional, however, in 
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that the existence of a liquid- vapour transition in ordinary 
dipolar system is still a hotly debated topic and one of the 
big unresolved questions regarding these particles [2H1I3S] • 
But as shown here, such a phase transition can be induced 
via an external driving field. 

The system considered in this study is driven and, con- 
sequently, inherently in a non-equilibrium state. The dy- 
namic coarsening observed in spinodal decomposition, on 
the other hand, is a process typically associated with non- 
driven systems. It is, however, not unique to those. For 
instance, active Brownian swimmers performing a "run- 
and-tumble" motion such as E. Coli bacteria, exhibit sim- 
ilar clustering behavior and demixing [5D]. With the on- 
going and rising interest in dynamics and non-equilibrium 
processes we expect an increasing amount of systems to be 
uncovered that are driven into cluster or pattern formation 
with behaviors similar to the one described here. 
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